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The Bardeen-Cooper-Schrieffer (BCS) theory for super- 
conductivity'' has been proved by several kinds of experi- 
ments. Among them, the isotope effect on superconducting 
transition temperature Tc has been one of key experiments 
for the BCS theory. If the Cooper pair is formed by phonon- 
mediated attractive interaction, Tc should be determined by a 
characteristic phonon energy cj, which is related to the mass 
of oscillator M as a; oc Af ^^/^. When we express the relation 
between Tc and M as 



Tr ex M-", 



(1) 



we obtain rj — 1/2 for BCS superconductors mediated by 
phonons. In actual experiments on Hg,^'^' it has been clearly 
shown that Tc is in proportion to M~^^^, leading to the ev- 
idence of phonon-mediated Cooper pair. Note, however, that 
in Ru,"*' the value of rj has been found to be smaller than 1/2. 
This is understood by the effect of Coulomb interaction in the 
famous McMillan formula.^' 

Recently, phonon-mediated superconductivity has attracted 
renewed attention from a viewpoint of anharmonicity since 
the discovery of superconductivity with relatively high Tc in 
/3-pyrochlore oxides. ^^' From both experimental and theoret- 
ical efforts, "^'"^' it has been gradually recognized that the su- 
perconductivity is induced by anharmonic oscillation of alkali 
atom contained in a cage composed of oxygen and osmium. 
Such anharmonic oscillation is called rattling and it is highly 
believed that the rattling plays a crucial role for the emergence 
of superconductivity in cage compounds. 

In this paper, we evaluate the exponent rj within the weak- 
coupling BCS theory for rattling-induced superconductor It 
is shown that 77 becomes larger than 1/2, indicating anoma- 
lous isotope effect. We also find that 77 increases with the in- 
crease of the amplitude of guest ion, which is relevant to /?- 
pyrochlore oxides. We propose an experiment on the isotope 
effect in order to prove a key role of rattling in /3-pyrochlore 
oxides. We emphasize that the increase of rj more than 1/2 is 
first reported in this paper, although the decrease of 77 less than 
1/2 has been understood by the effect of Coulomb interaction. 
Throughout this paper, we use such units as ft = fee = 1. 

We consider the Holstein model, given by 






H\oc 



(2) 



where k is momentum of electron, e^ denotes the energy of 
conduction electron, a is an electron spin, Cka is an annihi- 
lation operator of electron with k and a, H\oc denotes local 



electron-vibration term, expressed as 

H,oc = .g^ Qi4ci. + Y.^Pi/{'2M) + ViQi)]. (3) 

i,a i 

Here g is electron-vibration coupling constant, i denotes 
atomic site, Cia is an annihilation operator of electron at site 
i, Qi is normal coordinate of the oscillator. Pi indicates the 
corresponding canonical momentum, M is mass of the oscil- 
lator, and V denotes an anharmonic potential for the oscilla- 
tor, given by V{Qi) = Muj^Q'l/2 + k^Qj + kgQl Here w 
is energy of oscillator, while k^ and fcg are the coefficients for 
fourth- and sixth-order anharmonic terms. 

By using the phonon operator Ui defined through Qi = 
(at + ai)/V2AIuJ at site i, we obtain 



Hioc = \/a^y^(ai +at) + w^[atai + 1/2 

i,(T i 

+ P{ai + al)''+j{ai + al)% 



(4) 



where a = g^/{2APuj^), (3 = ki/{'iM'^uj^), and 7 = 
fc6/(8M'^a;^). With the use of non-dimensional parameters a, 
(3, and 7, it is convenient to rewrite V as 

auj{ql + l&aPqj + 640^75?), 



V{qi) 



(5) 



where qi is non-dimensional displacement, defined by qi = 
QiMuj^/g. 

Here we define the M dependence of parameters.'^' It is 
well known that the phonon energy lj is in proportion to 
j\//-i/2 from u! = y/k/M with a spring constant fc, when we 
assume that k dose not depend on M. If we further assume 
that g is independent of A/, we obtain a ex M^^^. Note that 
auj does not depend on A/. Concerning anharmonic param- 
eters /3 and 7, we obtain that /3 ex A/^^'^ and 7 ex A/^^ 
by assuming that fc4 and fcg are independent of M. Here- 
after, we explicitly consider ni dependence of parameters as 
Lo = ujo/y/m, a = ao\/m, /3 ~ fit^jy/m, and 7 = ^o/m, 
where m indicates the mass ratio of the guest ion. 

Now we consider the A/ dependence of the superconduct- 
ing transition temperature Tc, given in the BCS theory by 

1.13c.e-i/\ 



Tc 



(6) 



where A ~ U-p\^/W, C/ph is the phonon-mediated attraction, 
and W is the electron bandwidth. The exponent 77 in the iso- 
tope effect is evaluated by rj=—d log Tc/d log 771, leading to 

1 m dX 

Note that for a harmonic case with /3o=7o=0, we obtain 
X=2auj/W=2aouJo/W, which does not depend on A/. Thus, 
we find ?7 = 1/2 for the harmonic case (normal isotope ef- 
fect). However, for anharmonic phonons, [/ph depends on M 
and the value of 77 is deviated from 1/2. In the following cal- 
culations, we set W=l as an energy unit. 

First we show the anharmonic potentials considered in this 
paper. In Fig. 1(a), we show potentials for several values of 
/3o with 70=10""*, ajo=0.05, and Ao=2aowo=0.5, where j3'q = 
/3o/-v/7o- ^^ already mentioned in our previous paper,'"" the 
potential shapes are classified into three types: On-center type 
for i3'q > -V3/2, ratthng-type for -1 < ^^ < -\/3/2, and 



off-center type 



< -1. 
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Fig. 1. (Color online) (a) Anharmonic potentials for /3q=— 0.1, —0.866, 
— 1.0, and —1.1. Dotted curve denotes harmonic potential, (b) Variation 
of attraction (5f/ph/(^o ^s afunction of /9g for m=1.01. 
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Fig. 2. (Color online) (a) Exponent r] vs. /3' for 7=10^ , A=0.5, and 
id = 0.05. Dotted line indicates rj=l/2 for the normal isotope effect, 
(b) Exponent tj vs. Debye-Waller factor /. Note that the anharmonicity 
increases with the increase of /. 



Next we show that t/ph is actually changed by M for the 
case of anharmonic oscillation. For the purpose, we evaluate 
the magnitude of U-ph as 



t/, 



ph 
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e'^\ 
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where £'„ is the ground-state energy of H\oc for local elec- 
tron number n. In Fig. 1(b), we depict (5[/ph/wo as a function 
of /3q for 7o=10~^, Ao=0.5, and Ci;o=0.05, where ^f/ph is es- 
timated by (5f7ph = t/ph(m = 1.01) — Up\i{m = 1). For 
the diagonalization of iJioc, we use 250 phonon basis. Note 
that (5C/ph=0 for the harmonic case. We find that the attraction 
mediated by anharmonic phonons is significantly affected by 
M. '^^ In particular, for the on-center type potential, [/ph is in- 
creased by the increase of M, while it is decreased for rattling 
and off-center type potentials for 7o=10^"'. 

Now we move to the result of ry. In Fig. 2(a), we show 
?7 as a function of /3q for 7o=10~^, Ao=0.5, and ajo=0-05. 
Note that in the actual calculation, we evaluate dX/dm from 
SUpYi/O.Ol. For /3q > —0.8, we find the normal isotope effect 
with rj w 0.5, but in the region of —1.0 < /3q < —0.8 corre- 
sponding to the rattling-type potential, we find sharp increase 
of 7]. In the off-center type potential for (3q < —1.0, we find 
that 1] is moderately enhanced. 



Here we focus on the region of —1.0 < /3q < in or- 
der to discuss possible relevance of the present result with 
the isotope effect of /3-pyrochlore oxides AOS2O6 (A = K, 

Rb, and Cs). When alkali ion radius is decreased in the order 
of Cs, Rb, and K, the amplitude of the oscillation of alkali 
ion is increased, since the anharmonicity in the potential is 
increased in the order of Cs, Rb, and K."' Here we intro- 
duce the Debye-Waller factor / in a non-dimensional form 
as / = ((5?)/3£^, where £ == g/k and (•••) denotes the 
operation to take thermal average. In Fig. 2(b), we show ?] 
vs. /, where / is evaluated at T = 0.01, corresponding to 
a room temperature. Here we change anharmonic parameter 
set (/3, 7) from (-0.0253, 0.00451) for /3^=-0.35 and A=Cs 
to (-0.0392,0.00154) for /3o=-l-0 and A=K, which have 
been determined to reproduce the Debye-Waller factors for 
AOS2O6.''*' The anomalous isotope effect with 77 > 1/2 is 
expected to occur in /3-pyrochlore oxides. We predict that the 
value of rj is increased in the order of Cs, Rb, and K. 

In this paper, we have evaluated i] for rattling-induced su- 
perconductor on the basis of the BCS formula in the weak- 
coupling limit. As shown in our previous study, the strong- 
coupling effect is significant in the rattling-induced supercon- 
ductivity when the anharmonicity is increased. Thus, in order 
to discuss r; in more detail, it is necessary to calculate pre- 
cisely Tc as well as the renormalization factor by changing 
M within the Eliashberg theory. It is our future task. 

In summary, we have proposed that the isotope effect with 
the exponent 77 > 1/2 is found in the superconductivity due to 
electron-rattling interaction. The detect of this anomalous iso- 
tope effect can be the evidence of superconductivity induced 
by rattling in /3-pyrochlore oxides. 
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